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We investigate the impurity scattering rates for quasi-particles in vortex cores of sign-reversing 
s-wave superconductors. With use of the Born and Kramer-Pesch approximation for Andreev bound 
states, we show that the sign-reversed forward scatterings are dominant in vortex cores. Owing to 
the coherence factor in vortex cores of ±s-wave superconductors, the impurity scattering rate of the 
Andreev bound states has a characteristic distribution at the Fermi surfaces. For comparison, the 
impurity scattering rates in vortex cores of s-wave and d-wave superconductors are also discussed. 
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I. INTRODUCTION 

The discovery of novel Fe-based superconductors has 
attracted considerable attention because of high super- 
conducting transition temperature. 1 A ±s-wave pair- 
ing symmetry has been theoretically proposed as one 
of the candidates for the pairing symmetry in Fe-based 
superconductors The ±s wave symme- 
try means that the symmetry of pair potentials on each 
Fermi surface is s-wave and the relative phase between 
them is 7T. 

It is important for the identification of the ±s-wave 
symmetry to detect the sign change in the pair potentials 
between Fermi surfaces. Generally it is difficult to detect 
the relative phase of the order parameters in spatially 
uniform systems. As a phenomenon sensitive to the rel- 
ative phase, the formation of the Andreev bound states 
at a surface or junction has been discussed by several 
authors . 12 ' 13 The low energy excitations localized near 
vortex cores (Caroli-deGennes-Matricon modei^) can be 
regarded also as Andreev bound statesi 15 ' 16 ' 17 The im- 
purity scattering effect inside vortex cores is, therefore, 
expected to provide an efficient phase-sensitive tool. 

The effects of the phase-sensitive impurity scat- 
tering have been studied experimentally through the 
quasiparticle interference (QPI) by the scanning tun- 
neling microscopy/spectroscopy (STM/STS) in high- 
T c cupratesi 18 ' 19 i 20 : 21 i 22 : 23 For example, Hanaguri et 
al. have detected e?-wave coherence factors by the 
magnetic-field dependence of the QPI patterns in 
Ca2- a; Na 2 ;Cu02Cl2^ 22 With the use of this technique by 
the STM/STS, one might detect ±s-wave coherence fac- 
tors in the Fe-based superconductors^ 

In QPI experiments, the impurity scattering rates of 
quasiparticles in the absence and the presence of mag- 
netic field are necessary to identify the pairing symme- 
try. When the pairing symmetry in the Fe-based su- 
perconductors is ±s-wave, the energy spectrum is fully- 
gapped and hence the low energy excitations are ex- 
hausted by the Andreev bound states localized near vor- 
tex cores. One should consider the coherence effects of 



the Andreev bound states in the QPI patterns obtained 
by the STM/STS experiments. In earlier study^l how- 
ever, those effects have not been taken into account yet. 

The impurity scattering rates of the Andreev 
bound states in vortex cores have been treated 
analyticalf y 25 i 26 i 27 i 28 within the quasiclassical Eilcn- 
berger theor y 29 ' 30 for s-wav o 25 ' 26 and chiral p-wave 
superconductors ! 27 ' 28 Those analytical results on the im- 
purity scattering rates have been confirmed to be con- 
sistent with numerical results of self-consistent Born 
approximation ! 31 ' 32 The analytical calculation of the 
Andreev bound states in vortex cores originates from 
Kramer and Pesch^ in pure superconductors and has 
been generalized to the im pur e casei 25 ' 26 ' 27 We thus call 
the method used in Refs. I26II27I the Kramer-Pesch ap- 
proximation in this paper. 

The quasiclassical theory is applicable to Fe-based su- 
perconductor since the band widths (~ a few eV) are 
much larger than the maximum of the superconduct- 
ing gap (~ lOmeV). The iron-based superconductors are 
known to be multi-band systems and have multiple Fermi 
surfaces. 3 Assuming that the intraband pairings are dom- 
inant, one can define the quasiclassical Green functions 
on each band. One can then solve the equation of mo- 
tion of the quasiclassical Green function near a vortex 
core and can obtain the analytical expression of the im- 
purity scattering rate with the Born approximation in 
multi-band systems by extending the method developed 
in Refs. MM ■ 

The aim of this paper is to study impurity effects in 
vortex cores of two-dimensional multi-band superconduc- 
tors within the quasiclassical theory, with a particular 
interest in ±s-wave superconductors. We obtain the im- 
purity scattering rate from the initial state k to the 
final state k' by the Kramer-Pesch approximatio n 15 ' 26 ' 27 
in s-wave and d-wave and ±s-wave superconductors. 
We discuss the scattering-angle dependence of the sign- 
conserved and sign-reversed impurity scattering rates. 

This paper is organized as follows. The Born approx- 
imation and the quasiclassical approximation for multi- 
band superconductors are presented in Sec. II. The gen- 
eralization to multi-band superconductors is shown to 
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be straightforward when we assume that the intraband 
pairings are dominant. The impurity self-energy in a 
single vortex within the Kramer-Pesch approximation is 
derived in Sec. III. The analytical expression of the im- 
purity scattering rates m vortex cores is derived. 
The results are shown in Sec. IV. We discuss the coher- 
ence effects on the scattering-angle dependence for the 
cases of the sign-conserved and the sign-reversed scatter- 
ings. We also calculate the q = k' — fc-dependence of the 
impurity scattering rate for the isotropic s-wave, <i-wave 
and the isotropic ±s-wave superconductors, respectively. 
The conclusion is given in Sec. V. 



II. BORN APPROXIMATION AND 
QUASICLASSICAL APPROACH 

A. Orbital representation and Band representation 

Let us consider two-dimensional superconductors. We 
consider an n-orbital system which is a periodic crystal 
with n atomic orbitals in unit cell. We introduce the 
Hamiltonian written as 

H = eu c\ cu + At, ct s , +h.c. 



(1) 

where the operator c' k creates an electron with spin 
a and momentum k on the /x-th orbital. In the matrix 
form, this Hamiltonian is a 2n x 2n matrix in Nambu and 
orbital spaces expressed as 



H^(k x , k y 



H°(k x , ky) 

A°^k X ,ky) 



(2) 



in the "orbital representation" where the basis functions 
are atomic orbitals in crystal unit cell. Here, H° is the 
Hamiltonian in the normal state represented as annxn 
matrix in the orbital space and A° is the superconduct- 
ing order parameter. From now on, the subscript "o" 
indicates that matrices are represented with the orbital 
basis, hat a denotes an n x n matrix in the orbital space 
and check a denotes a 2n x 2n matrix composed of the 
2x2 Nambu space and the n x n orbital space. The 
unperturbed 2n x 2n Green function in the orbital rep- 
resentation is defined as 



(?o(fc x , ky\ mj„) = [iw n l - H^{k x ,k y )Y 



(3) 



Here, uj n is the Fermion Matsubara frequency. 

We also introduce an n x n Hamiltonian in the "band 
representation" defined by 

H h (k X ,ky) = P- 1 (k X ,k y )H (k X ,ky)P(k X ,k y ), (4) 



Ai 



(5) 



■• 

A„ 

Here, A; (i = 1, 2, • • • , n) denotes the i-th largest eigen- 
value. P is a unitary matrix consisting of the eigenvectors 



for the Hamiltonian H°(k x , k y ). The 2n x 2n Hamilto- 
nian in Nambu and orbital spaces in the "band represen- 
tation" is also defined by 



H^j(k x ,k y ) = U (k x ,ky)H^(k Xl ky)U(k x ,ky), (6) 



(7) 



where 



U(k x , ky ) 



P(k X ,ky) 



" P{k X ,ky) ' 



A b = P- 1 A°P. 



(9) 



The unperturbed 2n x 2n Green function in the band 
representation is written as 



Go{kxi ky] icj n ) — U(k x ,ky) GQ{k x ,ky\iui n )U{k x ,ky). 

(10) 

In general, A b contains off-diagonal elements, which cor- 
respond to inter-band pairings. 



B. Born approximation 

We consider the Green function under the influence 
of a lot of weak impurities. The Green function with 
the Born approximation in the orbital representation is 
written as 



G°(k) = G o (k) + G° (k)£°(k)G° (k), 



(11) 



where 
u{k) = 



dki 



(2nf 

u(k) 
-u(fc) 



u(k - fci)Gg(fci)fi(fci - fe),(12) 

(13) 



Here, ni mp is the concentration of impurity atoms, u(k) 
is the impurity potential. By substituting eq. (jTUJ) into 
eq. (fTTj) . the Green function in the band representation 
is written as 



G b (fc) = G b (k) + G h (k)£ h (k)G»(k), (14) 



where 



t b (k) = n imp / ^^fiCfc, fci)6S(fei)«(fci, fc), (15) 



(16) 



U " ! ' ; (27T) 2, 



v(k,ki) ee U(k)- x u(k-k x )U(ki) 



It should be noted that the impurity potential in the band 
representation v{k,k{) is a function of k and k\, not a 
function of k — k\. Assuming that intraband pairings 
are dominant, we neglect the off-diagonal (interband) el- 
ements in A b . In this case, the normal and anomalous 
parts of the self-energy E & ' N (fc) S b ' A (fc) are expressed as 
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dki 

/ (2^)2 



XI ^ m ( fe > fe O G 0,mm( fe l) V? fe), 

m 



(17) 
(18) 



where 



U(fe,fel) = P(fe) _1 M(fe - fcl)P(fcl), 



b _ 



Gb 




bt 



*8 

Gb 




(19) 
(20) 



Here, the symbols i,j,m denote the band indices. For 
simplicity, we assume that the perturbed Green function 
G h is diagonal in the band space. For example, the Green 
function in the two-band system is written a a 11 ' 33 ' 34 ' 35 



t G 



G = 



V 



F a \ 

Gp F p 

Fl G a 

-Fl Gp J 



(21) 



In this assumption, one can regard the self-energy as di- 
agonal with respect to the band index. 



J 



C. Quasiclassical Green functions 



We assume |Aj| <C E F . This relation is satisfied in 
most of systems such as conventional superconductors 
and Fe-based ones since the band width is a few eV 
and the superconducting order parameter is the order 
of lOmeV in these materials. In this case, one can use a 
quasiclassical approximation. 17 ! 29 ' 30 



Since the Green function Gq mm (ki) is localized 
around Fermi wave-length on the m-th band k™ F , the 
self-energy on the i-th band (fe) is written as 



s, b (fe) 



limp ^] 



dS 



F.k' 1 



J F,k" 



v im (k,ki) / d£ k ™G mm (ki) I v mi (ki,k) 



(22) 



r 



Here, k\ denotes the unit vector in the direction fci , v 
is the Fermi velocity on the m-th band and dS 



' F,k" 
F.k™ 1 



the Fermi-surface area element on the m-th band. We 
introduce the functions written as 



v(k F , k' F ) 

9{k F ) 
f(k F ) 



E S k F ,K A rn(kn (23) 
m 

X 5 k F ,k7 5 k' P ,kt imm ' ( fcf ' fc F ) ' ( 24 ) 

m,m' 

E*fe„fc? <f d^G b mm (k m ), (25) 
m J 

T, S k P ,k™ j <% km F* m (k m ), (26) 



where § d^k m shows that we take the contributions from 
poles close to the Fermi surface on the m-th band. With 
the use of the above functions, the normal and anomalous 



parts of the self energy are written, respectively, as 
£ N (fc F ) = n imp (v(k Fl k ltF )g(ki :F )v(k ltF ,k F )) FS , 



(27) 



£ (fc F ) = n imp (v(k F ,k^ F )f(k ltF )v(k hF ,k F ))FS, 



(28) 



(A(k' F )) FS ~ [ ^4 A (k' F ). (29) 



where, 



v F (k F ) 

Here we introduce the effective single-band Fermi surface 
as shown in Fig. [H since the relation fe™ k n F l is always 
satisfied for m ^ m' . Therefore, one can regard the ri- 
band system as a single-band system omitting the band 
index m and be allowed to translate a n x n matrix a 
into a scalar a from now on. 

In the quasiclassical approximation, the self-energy can 
be determined as a local value S(fep, r)M Therefore, one 
can calculate the self-energy by substituting g(k\,r) into 
Eq. (|27|) in inhomogeneous systems. 
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FIG. 1: Schematic figure of the effective Fermi surface. 



We introduce the quasiclassical Green function g de- 
fined by 



g{z,k Fl r) 



_( 9 f 



-f -9 J ' 



(30) 



which is a 2 x 2 matrix in the Nambu space and is a 
function of complex frequency z, Fermi wave-length k F , 
a point r = r(cos <j>, sin (f>) in real space. From now on, 
check a denotes a 2 x 2 matrix in the Nambu space. The 
equation of motion for g written as 

- iv F (k F ) ■ Vg = [zf 3 - A(r, k F ) -£,g], (31) 

supplemented by the normalization condition 

g 2 = -ir 2 L (32) 

Here, A is given by 



A(r,k F ) 



A{r,k F ) 



-A*(r,k F ) 







(33) 



and T3 is a Pauli matrix in Nambu space. In this paper, 
we consider the case where g is an analytic function of 
z in the upper half complex plane. Setting z = e + iS 
with real 5, we obtain the retarded Green function. We 
use a special parameterization of the quasiclassical Green 
function to solve the equation ^ The so- 

lution g of eq. pip can be written as 



«tt ( 1 — ab 2ia 

-2ib -(1-ab) 



1 + ab 



(34) 



Here, a and b are the solutions of the following Riccati 
differential equations: 

Of(fcjf)-Va = — 2(— iz + iXn)a 

- a 2 (A*-£ 21 ) + (A + £ 12 ), (35) 
v F (k F )-Vb = +2(-iz + ffiu)& 

+6 2 (A + £ 12 )-(A*-£ 21 ). (36) 

In this parameterization, the normalization condition 
(f3"2")l is automatically satisfied. For simplicity, we solve 
the Riccati differential equations under a given form of 
pair function. 



III. KRAMER-PESCH APPROXIMATION 

We calculate the self-energy in a single vortex core with 
the Kramer-Pesch approximation i 27 ' 42 i 43 Since Eqs. 



and (|36p contain V only through v F -~V, they reduce to a 
one-dimensional problem on a straight line, the direction 
of which is given by that of the Fermi velocity v F (k F ). 
We denote by a and b the unit vectors along the crystal 
axes X and Y, respectively. We introduce s, y, r by 



r = Xa + Yb, 

= si) + yu, 

r = Vx 2 + r 2 = Vs 2 + v 2 , 



with 



cos f„ sin t 
— sin 9„ cos ( 



(37) 
(38) 

(39) 



(40) 



Here, 9 V is the angle between a and the Fermi velocity 
v F (k F ). 

In the Kramer-Pesch approximation, we expand the 
coefficients of Riccati equations a and b with respect to 
the impact parameter y and the energy e. Near the vor- 
tex, the superconducting pair function in the first order 
with respect to y can be written as 

A(k F ,r) = /(IsDAooSignMfc^K 9 - (l + i-) .(41) 



A + Ai 



(42) 



Here, f(s) describes the spatial variation of the modulus 
of the pair potential and /(0) = 0, lim^oo f(s) = 1 
and d(k F ) describes the variation of the pair potential in 
momentum space. Following Refs. I27M42L we obtain the 
coefficients a and b: 



ai(s) 



with 



a = - sign(d(k F ))e l0 \ 
b = sign(d(k F ))e- w \ 

e «0) 



\v F (k F )\ 

e "( s ) 
\v F (k F )\ 



J da' (2ia (e - S) + 2i-^A 



ds 



(43) 
(44) 

,-u(«') 

(45) 

'(-2i6 (e-S)+2i|:At) e -«( s ') ) 

(46) 



E(s,fe F ) = Sii--sign(d(fc F ))(e^E 2 i+e-^S 12 ),(47) 



2\d(k F )\A x 

u(s) = — : — - — / ds tanh(s). 

\v F {k F )\ J 



(48) 



Here, we set f(s) = tanh(s). At a low energy e and a 
small impact parameter y, the Green function is written 
asSii 2 . 



-2m 



ai&o + aob_ 



-M, 



with 



M EE 



1 iaQ 
-ib -1 



(49) 



(50) 
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We obtain the impurity self-energy with the quasipar- 
ticle momentum k for the Kramer-Pesch approximation 
written as2£ 



-u(s) 



(51) 



since the denominator of g can be expressed as 



ai&o + ao»i = i._ jr. m I e - 2 y — ... m £fc 



|v F (fc F )| 



(52) 

Here, C = |« F (fc F )|/(2|d(fc F )|A oc ). 

Now we calculate the imaginary part of the impurity 



induced self-energy with the Born approximation. In 
pure superconductors (£& = 0), the Green function near 
a vortex core can be written as 

ir\v F (k F )\e-^ e-E(y,k F )-i6 - 
g (s,k F ,e) = — — ; — — : — ~ -M, 



C(k F ) (e-E(y,k F )) 2 + S 2 



where 



2e Ag o |d(fc F )| 2 

E[y, k F ) = y . 

\v F (k F )\ 



(53) 



(54) 



By substituting Eq. (|53|) into Eq. (|27|) . the self-energy 
expressed in Eqs. (l35|) and ([36]) becomes 



J 



Here, we introduce the coordinates (s",y') in the direction to the momentum fci >F : 

y = ssin(0„ - 6 v i) +y cos(0„ — 6 V >), (56) 
s" = scos(8 v -8 v <)-ysm(6 v -6 v >), (57) 

where 6 v i is the angle between a and the Fermi velocity u F (fei )F ). With the use of the above equation and 5 — > 0, 
ImS(s, fe F ) is written as 

Im£(s,y,fc F ) = n imv (v kF ki F V ki F kF - S(e-E(y', fei, F ))(l-sign[d(fe F )d(fei, F ))] cos(0„-<V))^ g . 

(58) 

Therefore, by performing the s'-integration in Eq. (|5ip . we finally obtain the imaginary part of the impurity induced 
self-energy: 

ImE fc = n imp (v k , k ,v k ,, k (1 - sign[d(fc)d(fc')] cos(0„ - M) ri ■ ff 2 n^^r^W 6 '"^^'^^)^. (59) 

Co| sin(6' w - 0„/)| |es(fc')| |«f(«)I 

with 

*<»■*> - ^^(™^-™k m(9 --^^ (60) 

s ' (k,k') = So (k, k') cos(e v - e v ,) ~ £ — M^L S in(e„ - g„>). (ei) 

The impurity scattering rate for the quasiparticles from the initial state k to the final state k' in a vortex core T ktk i 
is given by 

r M -(e) = n imp v k , k ,v k ,, k (1 - sign[d(fc)d(fc')] cos(0„ - 0„)) 71 • jf a ^SrW 6 ""' 80 ' 6 '^' ( 62 ) 

Co I sm(6»„ - |d(fc')| |«f(k)| 



r 



Then, the coherence factor for the impurity scattering in 
a vortex core is described by 

1 - sign[d(k)d(k')} cos(9 v ~ 9 V >). (63) 

The impurity scattering rate (|62p mainly depends on 
the impurity scattering matrix Vk,k'Vk',k^ the scattering 



angle (8 V — v i) and the sign-change of the scattering 
sign[d(k)d(k')]. 

The impurity scattering rate T k>k i also depends on the 
product of the local density of states for the quasiparticle 
with momentum k and that with k . The local density of 
states can be characterized by the impact parameters for 
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becomes zero: 



1 — cos(9 v — 6 v i) 
|sin(0„-0„,)| 



(65) 



In the case of the backward scattering (6 V — 9 V > = ir), 
in spite of the large coherence factor, the scattering rate 
is small since the point of the scattering is far from the 
vortex core so that the amplitude of the bound-state DOS 
becomes small as shown in Fig.^b). 




FIG. 2: Schematic figures of (a) the forward and (b) the back- 
ward scatterings around a single vortex core in the low energy. 



the quasiparticles. The quasiparticles with momentum k 
form a bound state around a vortex core. The impact pa- 
rameter y for these quasiparticles is defined from Eq. (|54")l 
and depends on the energy e and the magnitude of the 
Fermi velocity |t?f(fe)|. In the isotropic s-wave single- 
band superconductors, the pattern of the local density 
of states around a vortex core becomes circuler since the 
impact parameters for various quasiparticles are same42 
In the fully-gapped superconductors with multiple Fermi 
surfaces, there are several circuler patterns of the local 
density of states around a vortex core. At any case, one 
can clearly define the trajectory for the quasiparticles 
forming the Andreev bound states as shown in Fig. [2] 
Therefore, we can discuss the impurity scattering in vor- 
tex cores, in detail. 



B. Sign-reversed scattering 

In the case of d(k)d(k') < 0, the dependence of the 
scattering angle (9 V — 6 V >) in Eq. ([6"2")l can be written as 



1 + cos^ — 9 v t) 



sm 



(66) 



In contract with the sign-conserved scattering, the back- 
ward scattering (9 V — 9 v i = it) is suppressed because of 
the effect of the coherence factor. On the other hand, the 
forward scattering becomes quite large since the bound- 
state DOS becomes large at the scattering point near a 
vortex core as shown in Fig. [Ha). 

It should be noted that these sign-reversed forward 
scatterings hardly occur in single-band superconductors 
because quasiparticles close to each other in momentum 
space feel the same sign of the pair function. In multi- 
band superconductors with electron and hole circuler- 
like Fermi surfaces such as Fe-based materials, the inter- 
band scatterings can become the sign-reversed forward 
scatterings. The q-dependence of the impurity scatter- 
ing rate is anomalous in sign-reversing s-wave supercon- 
ductors since the intensity of the sign-reversed forward 
scatterings is much larger than that of all other scatter- 
ings. Here, q denotes the difference of the quasiparticlc's 
momentum q = k' — k. 



IV. RESULTS 



C. Numerical result 



A. Sign-conserved scattering 

We consider the case of d(k)d(k') > 0. In this case, the 
sign of the pair function does not change after the scat- 
tering. When the momentum dependence of the impurity 
scattering matrix Vk,k'Vk\k is not so large, the scattering 
angle (9 V — 6 V >) is important for the impurity scattering 
rate. The dependence of the scattering angle (9 V — 9 V >) 
in Eq. ([62]) can be written as 



1 — cos(9 v — v i) 



(64) 



One can find that the forward scattering (9 V — 6 V > = 
0) is completely suppressed since the coherence factor 



We show the q-dependence of the impurity scattering 
rate T^^+g in the two-band superconductors as a sim- 
plified model for Fe-based superconductors. We consider 
the system with two hole Fermi surfaces and two elec- 
tron Fermi surfaces as shown in Fig. [3] since Fe-based 
superconductors have the multiple Fermi surfaces^ We 
rotate crystal axes by 45°. The hole Fermi surface a\ 
(02) is located around (k x ,k y ) = (0,0) with the di- 
ameter 1.2/V2 (0.6/^/2) and the electron Fermi sur- 
face Pi (P2) is located around (k x ,k y ) — (tt/^/2, 7r/\/2) 
((k x ,k y ) = (-7r/V2,7r/\/2)) with the diameter 0.8/V2. 
For simplicity, the amplitude of the Fermi velocity is 
isotropic in momentum space |i)j?(fe)| = 1 and we set 
the impurity scattering matrix v^h^k' k by its average 
(vk.k' v k',k = !)• We consider various superconducting 
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FIG. 3: (Color online) Two hole Fermi surfaces ot\,a,2 and 
two electron Fermi surfaces /3i,/?2- The dashed (blue) square 
denotes the first Brillouin zone. 



(a) d-wave symmetry 




(b) ±s-wave symmetry 
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FIG. 4: (Color online) Schematic figures of the sign of the 
pair functions for (a):d-wave and (b):±s-wave symmetries. 




FIG. 5: (Color online) q-dependence of the impurity scatter- 
ing rate Yk,k+q in the isotropic s-wave superconductors. The 
energy is e = O.lAoo. 




FIG. 6: (Color online) q-dependence of the impurity scat- 
tering rate Y^.k+q m the d-wave superconductors d(k = 
2k x k v I '(kx + ky). The energy is e = O.lAoo. 



pairing symmetries; the isotropic s-wave d(k) = 1, the d- 
wave d(k) — 2k x k y /(k x + ky) and the isotropic ±s wave 
d(k) = — 1 for a bands and d(k) — 1 for /3 bands. 

First, we consider the isotropic s-wave superconductiv- 
ity. In this case, all scatterings are sign-conserved scat- 
terings. As shown in Fig. the impurity scattering rate 
is finite everywhere in q-space. This result is consistent 
with the result in Sec. IV. (A). 

Next, we consider the d-wave superconductivity as 
shown in Fig.[4ja). In this case, the scatterings between 
[3\ and /?2 Fermi surfaces are dominant since some of 
these scatterings are the sign-reversed forward scatter- 
ings. More specifically, the scattering from (k x ,k y ) = 
(— 7r + 0.8, 7r)/v2 (8 v = 0) on the fa Fermi surface to 
(k x ,k y ) = (-7T + 0.8, -tt)/V2 (9 v > = 0) on the fa Fermi 
surface is the sign-reversed forward scattering as shown 
in Fig.UJa). Here, we define the scatterings satisfying the 
relation v ~ v' as the forward scatterings. From Eq. (frlfi)) . 
one see that the impurity scattering rate diverges for this 
scattering. It should be noted that the condition for these 
sign-reversed forward scatterings to exist is sensitive to 
the shape of the Fermi surface and/or the direction of 
the Fermi velocity in the ri-wave superconductors. 

Finally, we consider the isotropic ±s-wave supercon- 



ductivity as shown in Fig. HJh). In this case, the inter- 
band scatterings are dominant. As shown in Fig. [3 
there are the arc-like strong intensity distributions in q- 
space. These intensity distributions are caused by the 
sign-reversed forward scatterings from the arc-like re- 
gions on the ot\ or«2 Fermi surfaces to those on the fa or 
fa Fermi surfaces. The regions where the sign-reversed 




FIG. 7: (Color online) q-dependence of the impurity scat- 
tering rate F k ,k+q in the isotropic ±s-wave superconductors. 
The energy is e — O.lAoo. 
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forward scatterings occur for the ±s-wave superconduc- 
tors are broader than those for the d-wave superconduc- 
tors. These properties are robust since the ±s-wave su- 
perconductors always have two or more Fermi surfaces. 
As shown in Eq. (|66p . one of the most important factors 
for the sign-reversed scatterings is the relation of the di- 
rections of the Fermi velocity of the quasiparticles before 
and after scatterings. In the systems with electron and 
hole Fermi surfaces such as Fe-based superconductors, 
the more the shapes of Fermi surfaces become similar to 
each other, the regions of the strong intensity become 
broader in q-space. 

One can easily extend our formulation into three- 
dimentional systems with use of the Fermi velocities pro- 
jected on the plane perpendicular to a vortex42 The Ric- 
cati equations (|35|) (|36|) for three-dimensional system turn 
into the same form as that for two-dimensional system 
because of a translational symmetry along the direction 
of the vortex. 

Our numerical results do not change qualitatively even 
when the gap-amplitudes on each Fermi surface are differ- 
ent or anisotropic. The impurity scattering rate depends 
on the anisotropy of gap-amplitudes only through the ra- 
tio of the gap-amplitudes between the initial states and 
final states |d(fc)|/|d(A;')| expressed in Eq. (|62|) . 

We discuss the possibility to observe the effects of the 
characteristic impurity scattering rate in vortex cores in 
±s-wave superconductors. The impurity scattering is 
one of the origins of the actual measurements by the 
STM/STS known as the QPI effects. It should be noted 
that the impurity scattering rate without vortex cores 
has not been understood clearly since the impurity scat- 
terings depend on not only the properties of supercon- 
ductivity but also the normal-states properties (i.e., the 
band-structure) and the type of the impurities. The 
STM/STS measurements can observe these effects as the 
QPI effects by applying magnetic fields. The results with 
vortex cores are more easily to understand. If the very 
strong arc-like peaks are observed and the intensity near 



q = is relatively small in the q-map of the dl/dV, 
this is one of the direct evidences of the ±s-wave super- 
conductivity. These peaks become large decreasing the 
energy. One should find the samples where the Andreev 
bound states can be observed more clearly so that the in- 
tensity of the sign-reversed forward scatterings becomes 
much stronger. 



V. CONCLUSION 

We studied the impurity effects in vortex cores of 
various kinds of superconducting pairing symmetries. 
We found that the sign-reversed forward scatterings are 
dominant for these impurity scatterings by the Andreev 
bound states in the low energy. The ±s-wave supercon- 
ductivity yields strong arc-like peaks and weak intensity 
near q = in the q-map of the dl/dV by the STM/STS 
measurements. We showed that it is important for the 
detection of the phase sensitive effects to investigate the 
system with magnetic fields. 
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